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I. INTRODUCTION 

The quark level flavor changing neutral interaction b — > sfi + fi~ is forbidden at 
the tree level in the standard model (SM) and can occur only at the one-loop level. 
Therefore it can serve as an important probe to test SM at loop level and also 
constrain many new physics models beyond the SM. This quark level interaction 
is responsible for the purely leptonic decay B s — ► and also the semi-leptonic 

decays B — > (K f K*)/j, + n~ . The semi-leptonic decays have been observed by BaBar 
and Belle [l|, y, u| with the following branching ratios: 

B(B -> Kfx + fT) = (5.7±l%) x 1(T 7 , 

B(B -> ICfi+fT) = (11.0ll;| 9 ) x iO -7 - (1) 

These values are close to the SM predictions BIB . However there is about 20% 
uncertainty in these predictions mainly due to the errors in the determination of the 
hadronic form factors and the CKM matrix element I V ts I ■ 



The decay B s — > is 
to be (3.35 ±0.32) x 10~ 9 



B J 

,7 



i^hly suppressed in SM. Its branching ratio is predicted 
, la, gI . This decay is yet to be observed experimentally. 
Recently the upper bound on its branching ratio has been improved to [h]] 

B(B S -> fi+fT) < 5.8 x 10" 8 (95% C.L.) , (2) 

which is still more than an order of magnitude above its SM prediction. B s — ► yU + /i _ 
will be one of the important rare B decays to be studied at the upcoming Large 
Hadron Collider (LHC) and we expect that the sensitivity of the level of the SM 



prediction can be reached with ~ 1 fb 1 of data. ll|, |12J. 



Many new physics models predict an order of magnitude enhancement or more in 
B(B S — ► These include theories with Z' mediated vector bosons |la^, as well 

as multi-Higgs doublet models that violate 13] or obey [ijj natural flavor conserva- 
tion. In [151 ]. it was shown that the new physics mediated by vector bosons is highly 
constrained by the measured values of the branching ratio of B — > (K, K*)fi + fi~. As 
a result, an order of magnitude enhancement in B(B S — > fi + from new physics 
vector or axial vector operators is ruled out. On the other hand, such an enhance- 
ment from the scalar /pseudoscalar new physics (SPNP) operators is still allowed, 



3 



since the most stringent bound on the SPNP operators comes from B(B S — > fi + 
itself. In particular, multi-Higgs doublet models or supersymmetric (SUSY) models 
with large tan/? can give rise to such an enhancement. 

Apart from the branching ratios of the purely leptonic and semi-leptonic decays, 
there are other observables which are sensitive to the SPNP contribution to b — > s 
transitions. These are forward-backward (FB) asymmetry App of muons [16) in 
B — > Kjx + n~ and longitudinal polarization (LP) asymmetry App of muons in B s — > 
[17]. Both these are predicted to be zero in the SM. Therefore, any nonzero 
measurement of one of these asymmetries is a signal for new physics. In addition, 
these asymmetries are almost independent of form factors and CKM matrix element 
uncertainties, which makes them attractive candidates in searches for new physics. 
In this paper we investigate what constraints the recently improved upper bound on 
B(B S — > ji^ //~) puts on the possible SPNP contribution to A FB and A LP . Do SPNP 
operators enhance these observables to sufficiently large values to be measurable in 
future experiments? 

The paper is organized as follows. In section [Til we study the effect of possible 
SPNP contribution to App. In section ITTl| we calculate the possible App enhance- 
ment due to SPNP, and point out some interesting experimental possibilities. In 
section [TV] we present our conclusions. 



II. FORWARD-BACKWARD ASYMMETRY IN B — > Kfi + ff 



There are numerous studies in literature of t 



SM and its possible extensions 18 



3, 20, E 



re FB asymmetry of leptons in the 



22 



23 



24). In the SM, the FB 



asymmetry of muons in B — ► Kfi + fi~ vanishes (or to be more precise, is negligibly 
small) because the hadronic current for B — > K transition does not have any axial 
vector contribution. However this asymmetry can be nonzero in multi-Higgs doublet 
models and supersymmetric models with large tan /3, due to the contributions from 
Higgs bosons. Therefore FB asymmetry in B — > Kfi + fi~ is expected to serve as an 
important probe to test the existence and importance of an extended Higgs sector 



21 



24j . Any nonzero measurement of this asymmetry will be a clear signal of new 



physics. 
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The average (or integrated) FB asymmetry of muons in B — > Kji^fi , which is 
denoted by (A FB ), has been measured by BaBar |2j and Belle |25|, |26| to be 



(A FB ) = (0.15±g;^ ± 0.08) (BaBar) , 



(3) 



(A FB ) = (0.10 ± 0.14 ± 0.01) (Belle). 



(4) 



These measurements are consistent with zero. But on the other hand, they can be 
as high as ~ 40% within 2a error bars. 



A. Calculation of A FB 

We consider new physics in the form of scalar /pseudoscalar operators. The ef- 
fective Lagrangian for the quark level transition b — > sfi + fi~~ can be written as 

L{b -> s[i + fi~) = L S m + L S p , (5) 



where 



(XG F \ rr 

JSM = ~l7r VtbV ts\ C 9 {si^Pib) fij^ + C 10 {S"f^P L b) //7 M 75/X 
- 2-^-m b (sia^Ppb) yq^ j , 

L SP = ^V tb V*\R s (sP R b)ilt, + R P (sP R b)p ll5 A . (7) 



V2ti 

Here Pl,r = (1 T 7s)/2 and is the sum of 4-momenta of /i + and fi~. Rs and Rp 
are the new physics scalar and pseudoscalar couplings respectively. In our analysis 
we assume that there are no additional CP phases apart from the single CKM phase. 
Under this assumption, Rs and Rp are real. Within SM, the Wilson coefficients in 
eq. have the following values: 

Cf = -0.310 , Cf = +4.138 + Y(q 2 ) , C 10 = -4.221 , (8) 



where the function Y(q 2 ) is given in 27J, [28]. 



The normalized FB asymmetry is defined as 

L dcosO , f T — f°, dcosQ , *f r a 

/ ^ JO azdcost) J —1 dzdcosv 

Jo dcos6 lzl£re + I-i dcosd 7 



dzd cos 9 
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In order to calculate the FB asymmetry, we first need to calculate the differential 
decay width. The decay amplitude for B{p) — > K{p')fi + {p + )fi~{p_) is given by 



x 



(K(p') \^b\ B(p)) {Cfu(p + )^v(p_) + C 10 u(p + )w 5 v(p_)} 



C cff 

-2-^m b (K(p>) \sia^q»b\ B(p)) u{p + )^v{p_) 

+ {K{p') \sb\ B{p)) {R s u(p + )v(p^) + R P u(p + ) l5 v(P-)} 



where — {p — p')^ = (p + + P-)^- The relevant matrix elements are 



, (10) 



(K(p>) \s^b\ B(p)) = (2p - q)J + (z) + (LJ!) qMz) - f + ( z )] , (11) 



(K{l/)\si<T fll ,</'b\B(i>)) =- (2p - q),q 2 - (m 2 B - m 2 K )q,\ h [ Z) , (12) 



{K{p') \sb\ B{p)) = mB{ \ e) f {z) . (13) 

17l b 

Here, k = mx/ms, z = q 2 /m 2 B and rhj, = m b /mB- In this paper, we approximate 
rh b by I. 

Using the above matrix elements, the double differential decay width can be 
calculated as 

^ ~ ° W \V th V;f m% ^ 2 {l,k 2 ,z)[3, 



dzdcosO 2 9 tt 5 



x 



|A| 2 ^ + | J B| 2 )^ + i0(l,A; 2 ,^)(|C| 2 + | J D| 2 )(l-^cos 2 
+2m M (l - k 2 + z)Re{BC*) + Am 2 |C| 2 
+2m M 0^(1, k 2 , z) 0^ Re{AD*) cos 9 



(14) 
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where 



.4 
B 



±(l-k 2 )f (z)R s , 



m^C 10 <j f+(z) - 
C = C w f + {z) 
D 

4_|_ J2 o, 7.2 /.2 

4m 



1 -^(fo(z) - f + (z))\ + i(l - k 2 )f (z)R P , 



C^'/ + (z) + 2C?"^, 
1 + /c 4 + z 2 - 2(fc 2 + /c 2 z + z) 

2 

1-^M . 



(15) 



Also, m M = m^/rriB and # is the angle between the momenta of K meson and n in 
the dilepton centre of mass frame. The kinematical variables are bounded as 

- 1 < cos# < 1 , 
4m 2 < z < (1 - k) 2 . 

The form factors /+,o,t can be calculated in the light cone QCD approach. Their 
q 2 dependence is given by 

f(z) = /(0) exp( Cl z + c 2 z 2 + c 3 z 3 ) , 



(16) 



where the parameters /(0),ci, c 2 and C3 for each form factor are given in Table [B 
The FB asymmetry arises from the cos 9 term in the last line of eq. (1140 . 





/(0) 


Cl 


C2 


C3 


u 


0.319lffi 


1.465 


0.372 


0.782 


fo 


o.3i9j5:8g 


0.633 


- 0.095 


0.591 


h 


n or:r:+0.016 
U.OOO_ 


1.478 


0.373 


0.700 



TABLE I: Form factors for the B ^ K transition 



18] 



The calculation of FB asymmetry gives 



A FB (z) 



2T m fl a 1 (z)(j ) (l,k 2 ,z)(3 2 R, 
dT/dz 



where 



r 2 nj 2 

u F u |2 5 



(17) 
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ai(z) = \{l ~ k 2 )C 9 f (z)f + (z) + (1 - k)C 7 f (z)f T (z) , (19) 

~ = \{l-k 2 )^zf Q {z){Rl + ^Rl) 

+ 2(1 - k 2 ) m M C 10 f (z)f + (z) ^ (z) (1 - k 2 + z)R P 

- 2(1 - k 2 ) m^Cxo^z^Mz) lf+(z) - l^(f ( z ) - /+(*))} R P 

+ 2mlC 2 10 P,4>kz) {/+(*) - ] -J^(f (z) - /+(*))' 2 

,2 ^2 



+ 8miCt (3^Hz)fi(z) 

1 2 777,^ 

+ _(i + —«i) W l( z )x 

O 2 



{A_r< e ff 2 Ar e 
(c?0+c ^ mw + i^_ /lw + !| +A) 

4m„ 2 C? /+(z) /3„ (1 - A; 2 + z) ^(z) x 



' 7 UWMz) 



h{z)- l -^{f {z)-f + {z))} . (20) 

From eq. ( 1TT1) . it is clear that Afb(z) is proportional to m M (~ 0.02), and to the 
scalar new physics coupling R$. In the minimal supersymmetric standard model 
(MSSM) and two Higgs doublet models, Rs itself is proportional to m M and tan 2 (3. 
Hence a large FB asymmetry is possible only for exceptionally large values of tan j3. 

The average FB asymmetry is obtained by integrating the numerator and de- 
nominator of eq. (|17p separately over dilepton invariant mass, which leads to 

_ 2T m fi pf l R s Jdza^z)^!,^^) _ 2r B T m M fi 2 R s J dza x {z) (p(l,k 2 ,z) 
{ FB} ~ T(B -> Kn+fj,-) ~ B(B -> Kfi+fi-) 

(21) 

where B(B — > Kjj^yr) is the total branching ratio of B — > Kfi + fi~ . The numerator 
in eq. ( 12TI) can be calculated to be 

2T B Y Q m tl 0f t R s J dza x (z)<\)(\,k 2 ,z) = (5.25 x 10~ 9 )(1 ± 0.20)i? 5 , (22) 

whereas the total branching ratio, including the contribution of SPNP operators, is 
given by [20| 

B(B -> Kfi + fi-) = [5.25 + 0.18(i2| + R 2 P ) - O.UR P ] (1 ± 0.20) x 10~ 7 . (23) 

In the SM calculation of BiB — > Kfi + two vector form factors, /o and / + , as 
well as the tensor form factor fa appear. The SPNP contribution, on the other hand, 
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G F 


= 1.166 x 1(T 5 GeV~ 2 


m Ba -- 


= 5.366 GeV 


a = 


7.297 x 1(T 3 


m B = 


5.279 GeV 


T B S 


= (1.437lffl) x 


Vtb = 


1.0 


T B d 


= 1.53 x 10" 12 s 


V ts = 


(40.6 ± 2.7) x 10~ 3 




= 0.105 GeV 


fBs = 


(0.259 ± 0.027) GeV [29] 


m K 


= 0.497 GeV 







TABLE II: Numerical inputs used in our 
taken from the Review of Particle Physics 



analysis. Unless explicitly specified, they are 



3C|. 



is only through /o- We have made the assumption that the fractional uncertainties 
in all the form factors are the same. The \V ts \ dependence in the numerator and 
denominator of eq. ( 12~TT) cancels completely, whereas the errors due to the form 
factors uncertainties cancel partially. We conservatively take the net error in (Afb) 
to be 30%, leading to 



(A FB ) 



5.25 x 10~ 9 R, 



[5.25 + 0.18(i2§ + R 2 P ) - 0.13R P ] x 10 



(1 ± 0.3) . 



(24) 



B. Constraints on (Afb) from B(B S —> ) 

We now want to see what constraints the present upper bound on B(B S — > fi + ^r) 
puts on the maximum possible value of (Afb)- The present experimental upper limit 
on B(B S — > /i + /U~) is an order of magnitude larger than the SM prediction. In such 
a situation, the SM amplitude for this decay will be much smaller than the new 
physics amplitude and hence can be neglected in determining the constraints on 
new physics couplings, Rs and Rp. In other words, we will assume that SPNP 
operators saturate the present upper limit. Therefore we need to consider only the 
contribution of L$p to the decay rate of B s — > fi + fi~. 
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The decay amplitude for B s — > fi + fi is given by 

olGf 1 
2y/2n 



s~1 

M(B a ^fi + fj,-) = ^^V tb V*{0\s^b\B s }[R s u(pMPii) + RpHP^Hv(Pii)] 



(25) 



On substituting 



we get 



(0 \s l5 b\ B s ) = -i iM^s- , (26) 



2V27T m b + m s 

(27) 

where m b and m s are the masses of bottom and strange quark, respectively. The 
calculation of the branching ratio B(B S — ► fi + fi") gives 

fi(fi s - /iV) = ^t^^ '^l 2 & + ^p) • (28) 

Here we have neglected terms of order m s /m b and approximated m_B s /mj by 1. 
Taking f Ba = (0.259 ± 0.027) GeV, we get 

B(B S -> = (1-43 ± 0.30) x 10~ 7 + R 2 P ) . (29) 

Equating the expression in eq. f[2"9l to the present 95% C.L. upper limit in eq. (TSJ), 
we get the inequality 

(R 2 S + R 2 P ) < 0.70 , (30) 

where we have taken the 2a lower bound for the coefficient in eq. (|29~I) . Thus, the 
allowed region in the Rs~Rp parameter space is the interior of a circle of radius 
~ 0.84 centered at the origin. 



In 3l|], it was shown that the SPNP operators cannot lower B(B —> Kfi + fi~) 
below its SM prediction. Therefore from eq. (124j) . the maximum value of (Afb) with 
the current upper bound on B(B — > Kfi + (jT) is 1.34% at 2a. If B(B S — > / u + /i~) is 
bounded to 10~ 8 , the 2a maximum value of (Afb) w iU De 0.56%. 

A naive estimation suggests that the measurement of an asymmetry (Afb) of a 
decay with the branching ratio B at na C.L. with only statistical errors require 
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number of events. For B — » Kfi + fi~, if (Afb) is 1% at 2cr C.L., then the required 
number of events will be as high as 10 18 ! Therefore it is very difficult to observe 
such a low value of FB asymmetry in experiments. Hence FB asymmetry of muons 
in B — > Kfi + fi~ will play no role in testing SPNP. 



III. LONGITUDINAL POLARIZATION ASYMMETRY IN B s -> fi + pT 

The longitudinal polarization asymmetry of muons in B s — > /i + yU~ is a clean 
observable that depends only on SPNP operators. It vanishes in the SM, whereas 
its value is nonzero if and only if the new physics contribution is in the form of scalar 
operator. Therefore any nonzero measurement of this observable Alp will confirm 
the existence of an extended Higgs sector. The observable Alp was introduced in 
ref . [3] , though the correspon ding analysis in the context of Kl —>■ had been 



carried out earlier 



32 



33 



34 



35]. In this section, we will determine the allowed 
values of Alp consistent with the present upper bound on B(B S — > fi + Li~), and 
explore the correlation between these two quantities. 

The most general model independent form of the effective Lagrangian for the 
quark level transition b — > sn + fi~ that contributes to the decay B s — > /U + /i~ has the 



form 



36 



37] 



+-Rs(s75&)GM) + Rp(s>y 5 b)(fi>y 5 iJ,)} , (32) 

where Rp, Rs and Ra are the strengths of the scalar, pseudoscalar and axial vector 
operators respectively. Note that the effective Lagrangian in eq. (|32l) is essentially 
the same as the effective Lagrangian given in eq. (jSJ). Here we have dropped C-j and 
Cg terms which do not contribute to B s — > fi + fi~ . In addition, the Ra in eq. (I32p 
is the sum of SM and new physics contributions. 

In SM, the scalar and pseudoscalar couplings i?| M and i?| M receive contributions 
from the penguin diagrams with physical and unphysical neutral scalar exchange 
and are highly suppressed: 

Rl M = i?| M oc ~ 1(T 5 . (33) 

m w 
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Also, -R^ = Y(x)/sin 6w, where Y(x) is the Inami-Lim function [38 ] 



Y(x) 



x 



x — 



+ 



3;r 



x — i ' (x — iy 

with x = (m t /M w f. Thus, Rf 1 ~ 4.3. 

The calculation of the branching ratio gives 



In a; 



(34) 
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36| 



B(B„ 



where 



m, 



m b + m s 



-Rv 



Am 2 :, 



m 



B s 



m 



B s 



m b + m s 



-R, 



64tt 3 



\Vt s Vtb | T B J 2 B m Bs Jl 



4ml 



m 



(35) 



(36) 



B a 



Here t# s is the lifetime of B s . Eq. fl35|) represents the most general expression for 
the branching ratio of B s — > fi + fi~. 

We now derive an expression for the lepton polarization. In the rest frame of 
we can define only one direction p the three momentum of \i~ . The unit 
longitudinal polarization 4-vectors along that direction are 



V 



(0, e 



0. 



V- 



V 



(37) 



Transformation of unit vectors from the rest frame of /i + to the center of mass frame 
of leptons (which is also the rest frame of B s meson) can be accomplished by the 
Lorentz boost. After the boost, we get 



V 



p_ 

TO, 



where is the muon energy. 

The longitudinal polarization asymmetry of muons in B s 

r(et) - r 



A 



LP 



net) + r(-e±) 



(38) 

is defined as 
(39) 



Thus we get 
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1 
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-0.6 




FIG. 1: A LP vs i? s plot for £(£ s 



(5.8,3.0,1.0) x 10 



-8 



with A\ p = A^ P = Aj^p. It is clear from eq. f!40p that Alp can be nonzero if and 
only if Rs ^ 0, i.e. for Alp to be nonzero, we must have contribution from SPNP 
operators. Within the SM, Rs — and hence Alp — 0. 

Using eq. (|35j) . we can eliminate Ra and Rp from eq. (jlOjl in favour of the physical 
observables B(B S — ► jU + /i~) and a s . We get [17] 



.4 



2a, 



LP 



4m? 



m 



m| sJ R 5 lB(B s ^fi+n-) 



4m? \ 
_ | 1 Ji ] 



m| sJ R 5 



m 6 + m s 



(41) 



m 6 + m s V a s V m z? s / 

Eq. (14 ip represents a general relation between the longitudinal polarization asym- 
metry Alp and the branching ratio of B s — > fi + fi~. 

We now explore the correlation between Alp and B(B S — > yU + /i~). It is quite 
obvious that when -B(-B s — > [i + fi~) ^ 10~ 8 , we can neglect the SM contribution in 
obtaining the bounds on i?s and -Rp. However if B(B S — ► fi + fi~) is of the order of 
the SM prediction, then we will have to take into account the SM contribution as 
well. Therefore it is reasonable to consider both the cases separately 
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1.0 2.0 3.0 4.0 5.0 6.0 7.0 

B(B s -> nV)x10 8 



FIG. 2: Plot between |^4lp| and B(B S — > pi + fi ) for different Rs values, when B{B S — > 
[i + u~) ;> 10~ 8 . The region B(B S — > [i + fi~) > 5.8 x 10 -8 is ruled out by experiments to 
95% C.L.. 

A. B(B S ^fi+fi-) % 10- 8 

We first consider the constraints on Alp coming from the present upper bound 
on B(B S — > Fig. [T] shows the plot between Alp and i?^ for three different 

values of B(B S — > ^ 10~ 8 . Fig. [2]is a plot between |v4^p| and -£>(£> s — > 

for various allowed values of The bands in Figs. [T] and [2] are mainly due to the 
uncertainties in CKM matrix element \Vts\ and decay constant Jb b - 

We see from Fig. [I] that the maximum possible value of Alp consistent with the 
present upper bound on B(B S — > is 100%, i.e. the present upper bound of 

B(B S — > does not put any constraint on Alp- Indeed, B(B S — > will 

be unable to put any constraint on Alp even if it is as low as 10~ 8 . 

Thus we see that the recently improved upper bound on the branching ratio of 
B s — > which provides the most stringent bound on SPNP couplings, fails 

to put any bound on Alp- Therefore Alp is more sensitive to SPNP operators as 
compared to B(B S — > fi + fi~). Any nonzero measurement of Alp will be evidence for 
an extended Higgs sector. 
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We would like to emphasize another important point: The measurement of 
B(B S — > /U + /i~) will only give the allowed range for the values of the SPNP couplings 
Rs and Rp. However the simultaneous determination of B(B S — > and Alp 

will allow the determination of new physics scalar coupling Rs (see Fig. [2]) and this 
in turn will enable us to determine the new physics pseudoscalar coupling Rp. 



B. B(B S 10" 



LHC is expected to reach the SM sensitivity in B s — ► fi + fj, . In fact, it may 



even go 5cr below the SM prediction ll[. Therefore it is worth considering the 



case when B(B S — ► yU + /i~) is of the order of the SM prediction. In this section we 
study the correlation between Alp and B(B S — > fi + fj>~) under the assumption that 
B(B S — > ji + fj,~) is close to its SM prediction. 
Taking R A = R S A M , eq. ([35]) gives 

B(B S -> p; + /0 = a s [28.8(i2| + #p) - 9.7i? P + 0.8ll , (42) 

which leads to 

R% + ( flp - 0.165)' = gjggBCgL^jO . (43) 

This corresponds to a circle in i?^ — Rp plane with centre at (Rs = 0, Rp = 0.165) 
and radius given by r = a/0.035 B(B s — > //+ /i~)/a s . 

Fig. [3] shows the plot between and i?s for three different values of -B(5 S — > 
/^ + A* - ) ^ 10~ 8 . Fig.HJis a plot between \ Alp\ and £>(-£?, — ► /i + /U~) for various allowed 
values of Rs- It is obvious from fig. [3] that Alp can be 100% even if B(B S — > /U + /i~) 
is close to its SM prediction. 

We now consider three exciting experimental possibilities, all of which can be 
accounted for with SPNP. 



1. B(B S —* fi + fx ) is consistent with SM but Alp 7^ 

It is possible to have a non-zero value of Alp even if B(B S — > fi + fi~) is equal to 
its SM prediction. We can re-write eq. (135j) in the following form: 

B(B S -> //+ /1-) = a s [(6 5M - &p) 2 + 6|] , (44) 
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FIG. 4: Plot between |Alp| and B(B S — ► fx + [i~) for different i?s values, when B{B S — ► 
[i + <^ 10~ 8 . The vertical shaded band corresponds to la theoretical prediction within 
the SM. The region B(B S — ► /i + //~) ^5x 10~ 9 is not ruled out; here we just concentrate 
on the region near the SM prediction. 
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where 
bsM 



rn 



B, 



-R f 



Am 2 , 



m 



B, 



m B rrib + m s 



-R-. 



rrib + m 

Here we have taken Ra = R& , 
the SPNP operators. Now if B(B S — ► /i + /i~) is equal to its SM prediction, then 



(45) 



Ra, i.e. we have considered new physics only through 



a s [(b S M ~ bp) 2 + b 2 s ] = a s b 



SM i 



which leads to 



(bi 



b S M? + bl 



h 2 

U SM ) 



or 



(m b + m s ) 

ftp 5 "SM 



m 



B, 



mi, + m s . 



m 



'SM 



(46) 



(47) 



(4J 



Eq. (|48|) represents a circle in Rs — Rp plane with center at (0, (to& + rn s )bsM/m 2 Bj- 
The circle representing eq. (1481) passes through the origin (Rg = Rp = 0), which 
corresponds to the SM. However, in general the points on the circle have nonzero 
Rs, and hence imply nonvanishing Alp- Therefore it is possible to have a nonzero 
value of A LP even if B(B S — > /U + /i~) is equal to its SM prediction. Thus A LP can 
still serve as an important observable to search for SPNP even if B(B S — ► /i + yU~) is 
observed to be very close to its SM prediction. 



2. LHCb fails to find B s — * /i + fi 

If LHCb fails to find B s —>■ fi + pT or puts an upper bound on its branching ratio 
which is smaller than 2 x 10~ 9 (5a below SM prediction), this scenario can still be 
accomodated within the SPNP. 

The interference between the SPNP and SM operators can decrease the branching 
ratio B(B S — > /i + /i~) far below its SM prediction. In fact it can be seen from 
eq. fl35|) . B(B S — ► /U + /i~) can even vanish, provided the following conditions are 
satisfied simultaneously: 

R s = 0, R P = 2m/ f" b fl A = 0.04 R A . (49) 

m B s 

From Fig. HI it can be seen that for low Rs values, it is indeed possible to suppress 
B(B S — > ji + /i~) much below its SM value. 
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3. Both B(B S — ► /i + ) and Alp are consistent with the SM 

The lepton polarization asymmetry is a result of the interference of the scalar 
term with pseudoscalar / axial vector, as can be seen from eq. (HU1) . Therefore it 
vanishes when either bs or (bp — &sm), defined in eq. (H5l) . vanishes. Thus there exists 
the interesting possibility of nontrivial SPNP even when both B(B S — > [i + and 
Alp are consistent with the SM. This occurs when 

b s = 0, b P = 2b SM , (50) 
or b s = b S M , bp = b SM , (51) 

as can be confirmed from eq. (H71) . Therefore, the absence of SPNP is not guaranteed 
simply by the consistency of these observables with the SM; more channels need to 
be examined to rule out this possibility completely. 

IV. CONCLUSIONS 

An order of magnitude enhancement in B(B S — > /U + /i~) is possible only due to 
SPNP operators. Apart from B(B S — > /i + yU~), observables such as FB asymmetry 
of muons in B — > Kfi + \l~ and LP asymmetry of muons in B s — > are also 

sensitive to SPNP operators. In this paper we consider the constraints on possible 
SPNP contribution to these observables coming from the present upper bound on 
B(B S -^fi+fT). 

We find that B(B S — > fi + puts very stringent constraint on SPNP contribution 
to (Afb) an d restricts its value to be less than ~ 1%. Such a small FB asymmetry 
is almost impossible to be measured in experiments. In the literature, (Afb) of 
muons in B — > Kfi + fi~ has been considered a promising measurement for probing 
SPNP operators. Our results show that the present upper bound on B(B S — > ji^) 
makes searching for SPNP through (Afb) a futile exercise. 

On the other hand, the present upper bound on B(B S — > fi + Li~) does not put 
any constraint on A L p. Indeed, A L p can be 100% even if B(B S — > /i + /U~) is close to 
its SM prediction. Alp is sensitive only to SPNP operators and hence its nonzero 
value will give direct evidence for a non-standard Higgs sector. 
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A simultaneous determination of B(B S — > ^ + ^) and will enable us to 
separate the new physics scalar and pseudoscalar contributions. Therefore it is worth 
considering this observable in experiments to probe the extended Higgs sector. 
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